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Abstract 



O 

^■i We report a concise answer-in the case of 2 x 2 systems-to the fundamental quantum-information- 



theoretic question as to "the volume of separable states" posed by Zyczkowski, Horodecki, Sanpera 
and Lewenstein (Phys. Rev. A, 58, 883 [1998]). We proceed by applying the Mathematica 
command FindSequenceFunction to a series of conjectured Hilbert-Schmidt generic 2x2 (rational- 
valued) separability probabilities vr(a), a = 1,2,... ,32, with a = 1 indexing standard two-qubit 
systems, and a = 2, two-quater(nionic)bit systems. These 32 inputted values of 7r(a)-as well as 
32 companion non-inputted values for the half- integers, a = \ (two-re [al] bit systems), |,...^- 
are advanced on the basis of high-precision probability-distribution-reconstruction computations, 
employing 7,501 determinantal moments of partially transposed 4x4 density matrices. The func- 
tion P(a) given by application of the command fully reproduces both of these 32-length sequences, 
and an equivalent outcome is obtained if the half-integral series is the one inputted. The lengthy 
expression (containing six hypergeometric functions) obtained for P{ct) is, then, impressively con- 
densed (by Qing-Hu Hou and colleagues), using Zeilberger's algorithm. For generic (9-dimensional) 
two-rebit systems, P{\) = §§, (15-dimensional) two-qubit systems, P(l) = (27-dimensional) 

& '. 

two-quaterbit systems, P{2) = while for generic classical (3-dimensional) systems, -P(O) = 1. 
PACS numbers: Valid PACS 03.67.Mn, 02.30. Zz, 02.30. Gp, 02.40.Ky 
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FIG. 1: Generalized two-qubit separability probability function P(a), with P(0) = l,P(\) = 
||,P(1) = J|,P(2) = ^ for generic classical four-level (a = 0), two-rebit (a = two-qubit 
(a = 1) and two-quaterbit (a = 2) systems, respectively. 

We present a succinct formula P(a) of fundamental quantum-information-theoretic im- 
portance. It effectively addresses, as well as considerably expands upon, the 2x2 instance of 
the question posed by Zyczkowski, Horodecki, Sanpera and Lewenstein in their highly-cited 
1998 paper, "Volume of the Set of Separable States" These authors gave "three main 
reasons of importance" -philosophical, practical and physical — for studying the question of 
"how many entangled or, respectively, separable states there are in the set of all quantum 
states". The concise formula P{pi) takes the form (Fig. [[]) 

P(a) = £~ /(a + z), (1) 

where 

/(a) = P(a) - P(a + 1) 

and 

q(a) = 185000a 5 + 779750a 4 + 1289125a 3 + 1042015a 2 + 410694a + 63000 = (3) 
a(5a(25a(2a(740a + 3119) + 10313) + 208403) + 410694) + 63000. 

Now, we claim that the specific values of interest given by P{\) = ||,P(1) = J,, and 
P(2) = ^ are the separability probabilities (with respect to Hilbert-Schmidt measure 




g(a)2~ 4a - 6 r(3a + §)r(5a + 2) 
3r(a + l)r(2a + 3)r(5a + f ) 



(2) 
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of generic (9-dimensional) two-rebit systems, (15-dimensional) two-qubit, and (27- 
dimensional) two-quater(nionic)bit systems, respectively (cf. fl). Further, P(0) = 1, for 
generic (3-dimensional) classical systems. 

The derivation of these results-computationally/numerically strongly convincing, though 
still not formally rigorous in nature-was begun in the paper "Moment-Based Evidence for 
Simple Rational- Valued Hilbert-Schmidt Generic 2x2 Separability Probabilities" ^J. Two 
basic sets of analyses were conducted there. 

The first set was concerned with establishing formulas for the bivariate determinantal 
product moments ( |p PT | n |p| / ,k,n = 0, 1, 2, 3, ... , with respect to Hilbert-Schmidt (Eu- 



clidean/flat) measure {2] {3, sec. 14.3], of generic two-rebit (4, |6[ and two-qubit density 
matrices (p). Here p PT denotes the partial transpose of the 4x4 density matrix p. Nonneg- 
ativity of the determinant |p PT | is both a necessary and sufficient condition for separability 
in this 2x2 setting [7J. 

In the second set of analyses in [5J, the univariate determinantal moments (|p PT |") 
and ((|p PT | |p|) n ), induced using the bivariate formulas, served as input to a Legendre- 
polynomial-based probability distribution reconstruction algorithm of Provost [a, sec. 2]. 
This yielded estimates of the desired separability probabilities. (The reconstructed proba- 
bility distributions based on |p PT | are defined over the interval |p PT | G [— ^, jjlg], while the 
associated separability probabilities are the cumulative probabilities of these distributions 
over the nonnegative subinterval |p PT | G [0, ^g]. We note that for the fully mixed (classical) 
state, |p PT | = 2I3, while for a maximally entangled state, such as a Bell state, |p PT | = — A.) 

A highly-intriguing aspect of the (not yet rigorously established) determinantal moment 
formulas obtained (by C. Dunkl) in [5, App.D.4] was that both the two-rebit and two-qubit 
cases could be encompassed by a single formula, with a Dyson- index- like parameter a {9] 
serving to distinguish the two cases. The value a = \ corresponded to the two-rebit case 
and a = 1 to the two-qubit case. Further, the results of the formula for a = 2 and n = 1 and 
2 have recently been confirmed computationally by Dunkl using the "Moore determinant" 
(quasideterminant) 10|, 111] of partially transposed 4x4 quaternionic density matrices. 

When the probability-distribution-reconstruction algorithm jg] was applied in j^] to the 
two-rebit case (a = ~), employing the first 3,310 moments of |p PT |, a (lower-bound) estimate 



2/> 

that was 0.999955 times as large as || ~ 0.453120 was obtained (cf. [12j, bot. p. 6]) 



Analogously, in the two-qubit case (a = 1), using 2,415 moments, an estimate that was 
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0.999997066 times as large as J| m 0.242424 was derived. This constitutes an appealingly 
simple rational value that had previously been conjectured in a quite different (non-moment- 
based) form of analysis, in which "separability functions" had been the main tool employed 

m. 

Further, the determinantal moment formulas advanced in jsj were then implemented 
with a set equal to 2. This appears to correspond-as the indicated recent computations of 



15|- 



Dunkl support-to the generic 2 7- dimensional set of quaternionic density matrices 
Quite remarkably, a separability probability estimate, based on 2,325 moments, that was 
0.999999987 times as large as g « 0.0804954 was found. 

In this study, we extend these three (independently conducted) moment-based analy- 
ses in a more systematic manner, jointly embracing the 64 integral and half-integral values 
a — |, 1, §, 2, . . . , 32. We do this by accelerating, for our specific purposes, the Mathematica 
probability-distribution-reconstruction program of Provost jg], in a number of ways. Most 
significantly, we make use of the three-term recurrence relations for the Legendre poly- 
nomials. Doing so obviates the need to compute each successive higher-degree Legendre 
polynomial ab initio. 

In this manner, we were able to obtain-using exact computer arithmetic throughout- 
" generalized" separability probability estimates based on 7,501 moments for a = 
|, 1, |, . . . , 32. Now, the associated two-rebit separability probability estimate in- 
creased/improved from 0.999955 times as large as || « 0.453120 to 0.999989567; the 
two-qubit separability probability estimate improved from 0.999997066 times as large as 

~ 0.242424 to 0.99999986; and the two-quaterbit separability probability estimate im- 
proved from 0.999999987 times as large as J§ « 0.0804954 to 0.999999999936. 

We note that the zeroth-order approximation (being independent of the particular value 
of a) provided by the Provost probability-distribution-reconstruction algorithm is simply the 
uniform distribution over the interval [—jq, ^k]- The corresponding zeroth-order separability 
probability estimate is the cumulative probability of this distribution over the nonnegative 
subinterval [O,^], that is, ^/(ici 2^) = Tf ~ 0.0588235. So, it certainly appears 
that speedier convergence of the algorithm occurs for separability probabilities, the true 
values of which are initially close to ^ (such as ^ ~ 0.0804954 in the quaternionic case). 
Convergence also markedly increases as a increases. 

It appeared, numerically, that the generalized separability probabilities obtained for all 
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the 32 integral and 32 half-integral values of a employed were rational values (not only for 
the three specific values a — |, 1, 2 of original focus). With various computational tools and 
search strategies based upon emerging mathematical properties, we were able to advance 
additional, seemingly plausible conjectures as to the exact values for the remaining 61 values 
of a, as well [la, p. 6]. 

We fed the sequence of 32 integral-indexed conjectured rational numbers into the Find- 
SequenceFunction command of Mathematica. ("FindSequenceFunction[list] uses earlier ele- 
ments in list to find candidate simple functions, then validates the functions by looking at 
later elements. FindSequenceFunctionflist] only returns functions that correctly reproduce 
all elements of list.") This produced a generating formula P(u) incorporating a diversity of 
hypergeometric functions of the p F p _i type, p = 7, . . . , 11, all with ar gum ent z = |j = (|) 3 . 



(We note that z a = <i/§t is the "residual entropy for square ice" 



13, P- 412] (cf. 
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eqs.[(27), (28)]) [19].) In fact, the Mathematica command succeeds using only the first 
twenty-eight conjectured rational numbers, but no fewer-so it seems fortunate, our compu- 
tations were as extensive as they were. The formula P(a) produced was quite cumbersome 



in nature (extending over several pages of output) (cf. |16|, Fig. 3]). 



P(a) did yield values for half- integral a (including the a — % two-rebit conjecture of ||), 
also exactly fitting our corresponding 32 half-integral rational- valued conjectured separabil- 
ity probabilities [16J, p. 6]. The process was fully reversible, with the half- integral sequence 
serving as the inputted one. 

We had been unable to find an equivalent form of P(ct) with fewer than six hypergeometric 



functions 
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Fig. 3]. Qing-Hu Hou and colleagues, however, were able to obtain the 
remarkably succinct results ([I])-®. They did so, first observing that the hypergeometric- 
based formula for P(a) could be expressed as an infinite summation. Letting Pi(a) be 



the l-th such summand, application of Zeilberger's algorithm 20| (a method for producing 
combinatorial identities) (cf. J21} App. A]), then, yielded that 

Pi{a) - Pi(a + 1) = -Pi +l {a) + P^a). (4) 



(The package APCI-available at http://www.combinatorics.net.cn/homepage/hou/-was 



employed.) Then, summing over I from to 00, they found that 

P(a)-P(a + l) = P (a). (5) 
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Letting f(a) = Po(a), the concise summation formula ([I]) is obtained. (C. Krattenthaler 
indicated-and Hou agreed-that these results might equally well be derived without recourse 
to Zeilberger's algorithm.) 

We certainly need to indicate, however, that if we do explicitly perform the infinite sum- 
mation indicated in (pQ), then we revert to a ("nonconcise") form of -P(a), again containing 
six hypergeometric functions. Further, it appears that we can only evaluate ([1]) numerically- 
but then easily to hundreds and even thousands of digits of precision-giving us extremely 
high confidence in the specific rational- valued Hilbert-Schmidt separability probabilities ad- 
vanced. 

There remain the important problems of formally verifying the formulas for P{a) (as 
well as the underlying determinantal moment formulas in [5|, employed in the probability- 
distribution reconstruction process), and achieving a better understanding of what these 
results convey regarding the geometry of quantum states [3]. Further, questions of the 
asymptotic behavior of the formula (a — > oo) and of possible Bures metric 3], 
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terparts to it, are under investigation [26]. 



25l | coun- 



The foundational paper of Zyczkowski, Horodecki, Sanpera and Lewenstein," Volume of 
the set of separable states" [l|], did ask for volumes, not specifically probabilities. At least, 
for the two-rebit, two-qubit and two-quaterbit cases, a — |, 1 and 2, we can readily convert 
the corresponding separability probabilities to the separable volumes g^fff^ = 2^I^W' 

_ QnH ill — Till lisinp- thp 

449513064000 2 6 -3 6 -5 3 -7 2 ll 2 13 3914156000371803494400000 2 14 ?. ln -5 5 7 3 ll 2 13 17 2 -19 2 -23 ' UD111 6 



Hilbert-Schmidt volume formulas of [14, Thms. 1-3] (cf. [2|, yj). The determination of 
separable volumes-as opposed to probabilities-for other values of a than these three appears 
to be rather problematical, however. 



Let us also note that Theorem 2 of 



271 ] . in conjunction with the results here, al- 



lows us to immediately obtain the separability probabilities of the generic minimally- 
degenerate/boundary 8-, 14-, and 2 6- dimensional two-rebit, two-qubit, and two-quaterbit 
states, as one-half (that is, and ^) the separability probabilities of their generic 

non-degenerate counterparts. 

In addition to P(a) assuming rational values for nonnegative integral and half-integral 
values of a, we note that P(-\) = = 2,P(-1) = f and P(-|) = P{-\) = §. 



A number of other interesting (irrational) values of P(a) are indicated in |16|, p. 9]. For 
■ _ , - 3r £.y , where the term 



example, P{—\) = 2-\ — . 2 > where the term f \' « 1.46099848 is "Baxter's four-coloring 
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constant" for a triangular lattice 17], p. 413]. The first derivative P'(a)-similarly to P{oi) 



itself-appears to assume rational values for nonnegative integral and half-integral values of 

130577 
' 457380 ' 



a, such as P'(l) = -|pgg. We further observe that P'(0) = -2 and P"(0) = 40 -20C(2) 



40 — ~ 7.10132. Thus, let us conclude that the function P(a) certainly appears to be 
an object well worthy of further investigation. 
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